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Abstract:  Using  a log  normal  model  for  molecular  weight  distribution, 

drag  curves  can  be  predicted  which  agree  with  experiments  in 
pipe  flow.  Working  back  from  pipe  flow  data,  average 
molecular  weights  or  concentrations  can  be  estimated  from 
drag  curves.  Illustrations  are  given  for  a recent  v.’ater 
tunnel  tf;st  and  data  in  the  literature. 
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THE  INTERPRETATION  OF  POLYOMETER  DATA 


Introduction 


Drag  reduction  by  dilute  polymer  solutions  is  easily  Illustrated 
by  pressure  drop  vs.  flow  rate  curves  for  turbulent  flow  In  a circular 
pipe.  It  Is  desired  to  compare  such  tests  to  experiments  In  other 
geometries  with  polymers  of  different  molecular  weight  distributions 
and  different  solution  concentrations.  For  example  In  a polymer  ocean 
test  the  polymer  Initially  is  fresh  with  a certain  molecular  weight 
distribution,  and  as  it  degrades  it  may  be  equivalent  (for  drag 
reduction  at  least)  to  a fresh  solution  of  a lower  concentration. 

In  tl»e  analysis  In  this  appendix  we  first  describe  the  presentation  of 
pipe  flow  data  and  then  develop  a model  to  relate  the  curves  to 
concentration  and  molecular  weight  distribution. 


Presentation  of  Data 

Over  the  range  of  Reynolds  numbers  lO’’  to  10^  the  flow  rate 
for  turbulent  flow  can  be  related  to  the  shear  rate  by  a power  law 
relationship.  In  terms  of  Q In  cubic  cenclmetcrs  per  second  and 
U^^/v  In  sec**^  i 


0,572 
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Here  D is  the  pipe  diameter  and  U*/v  equals  ?w/p  equals 
(-Ap)D/4Lu  where  v is  the  kinematic  viscosity,  p,  the  viscosity, 
the  shear  stress  at  the  wall,  and  -AP/L,  the  pressure  drop 
per  unit  length.  The  limiting  drag  rcductlt ...  asymptote  as  given 
by  Virk  (1967)  can  also  be  expressed  in  n similar  form: 


- 0.71  d2 


.31 


Now  when  a curve  of  Q vs.  U*^/v  is  measured  using  the  polyometer 
described  in  another  section,  several  different  types  of  behavior 
are  found.  These  are  Illustrated  on  fig.  1.  Curve  I shows  a 
polymer  solution  containing  sufficient  polymers  of  a high  enough 
molecular  weight  to  exhibit  maximum  drag  reduction  at  all  P.eynolds 
numbers.  When  the  polymer  molecular  weight  is  the  name  as  for 
Curve  I but  the  concentration  is  too  low  to  reach  the  maximum 
drag  reduction  line.  Curve  II  results.  Keep  in  mind  that  the 
polymer  solutions  i>e  are  discussing  have  a molecular  weight 
distribution  and  are  not  of  a single  molecular  weight.-  The  next 
curve.  III,  illustrates  a solution  with  lower  molecular  weight 
polymers  than  in  Curves  I and  11  but  enough  concentration  to  reach 
maximum  drag  reduction  at  a sufficiently  high  value  of  shear  rate. 
In  this  case  the  curve  approximately  follows  the  vrnter  turbulent 
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line  up  Co  an  "onset  shear  sCrcas  and  Chen  proceeds  at  a f,rcaCer 
slope  to  Che  maxlinum  drag  reduction  line.  In  reality  the  deviation 
from  the  turbulent  water  and  the  intersection  with  the  tnaxlmum 
drag  reduction  lines  are  not  sharp  as  shown  but  very  gradual.  In 
almost  all  cases,  however,  the  curve  can  be  well  approximated  by 
a straight  line  drawn  between  the  water  and  maximum  drag 
reduction  lines.  Curve  IV  shows  another  such  line  with  not  enough 
experimental  data  at  very  high  shear  rates.  This  curve  would 
eventually  reach  the  maximum  drag  reduction  line.  Finally  Curve  V 
represents  the  curve  for  the  minimum  concentration  of  polymer  so 
Chat  at  a shear  rate  of  infinity  the  maximum  drag  reduction  curve 
Is  reached.  Other  curves  beyond  this  would  never  reach  maximum  drag 
reduction. 


Description  of  Data  ' 

Now  to  describe  these  curves  mathematically  we  can  consider 
only  those  of  types  III,  IV,  and  V.  Curve  I may  of  course  change  to 
Curve  III  if  a larger  size  pipe  is  used  in  the  polyoraeter  (lower 
transition  shear  stress).  We  can  define  two  types  of  fractions  to 
represent  the  points  between  the  turbulent  water  and  the  maximum 
drag  reduction  curves  as  follows: 


H Qp  **  Qw 


. ^ U*^/v  ■ Const. 
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I 

• » I 

fi 


fw  - fp 


■ Const. 


where  the  subscript  p Is  for  polymer,  w,  for  water,  and  M,  for 
maximum  drag  reduction,  and  f is  the  friction  factor.  For  the 
polyometer  curves 


(^)w- 

W).  - 


“ m. 


Const. 


Figure  2 shows  some  possible  and  Aj^  locations.  About  point  P 


we  have 


where  XY  denotes  the  number  X minus  the  munher  Y and  not  the  length 
or  log  coordinates.  For  many  x«tcr  tunnel  tests  we  need  a fractional 
drag  reduction  to  aid  in  the  explanation  of  the  data.  For  skin 
friction  drag  in  turbulent  flow,  when  the  water  would  be  at  point  B, 
the  fractional  drag  reduction  for  the  Illustrated  polymer  solution 
at  the  same  free  streem  velocity  is  (1-G/B).  If  an  experiment 
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involves  an  Irrotatlonal  flow  and  the  strain  rate  Is  the  Important 

I 

parameter  then  the  appropriate  functional  drag  reduction  for  the  j 
illustrated  polymer  for  a shear  rate  on  the  line  APB  is  (1-P/D). 


Log  Normal  Model 

We  would  like  to  relate  the  point  P and  all  points  on  the 
line  OPT  to  the  molecular  weight  distribution,  concentration,  and 
behavior  of  polymer  in  the  turbulent  flow.  Although  the  molecular 
weight  distribution  in  the  dilute  polymer  blend  is  irobably 
complex,  we  choose  a log  normal  distribution  as  representative. 
Then 


_ (inM  - ro)^ 


P(M)  - _1_  e 
/ZnBM 


j o2 


(6) 


where 

m ■ <fnM> 


and 


<M> 


<x> 


m + JsB2 

■ e 


" / 


o /TiTeM 


dM 


(7) 
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in  units  of  decilitere/grara.  Then  at  25®C  for  water  aa  the 
solvent 


- 2 X 10‘*5  j(1.78 


(12) 


To  relate  to  the  polyometer  curve  we  need  the  relationship  between 
Tj  and  the  flow  timescale  U*2/v.  Lumley  finds  that  the  condition 
for  ultimate  expansion  o;  the  polymer  molecules  in  turbulent  flow 
in  given  by 

2|xP'T>+j  (13) 


where  X is  related  to  the  shear  raid  for  the  pipe  flow  and 

is  an  Integral  time  scale.  Using  Lumley' a definitions  of  X and  ‘TT’ 

for  pipe  flow  we  obtain, 


Wlien  eq.  9 is  put  in  terras  of  T^^, 


(14) 


5S 
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6p  - 6w 


, -(JtnT,*5fi2  ^ -Count.) 2 

j“e  26^^ 

3 (Con8t)T0-5G2 


£Ii 

Tj.56? 


(15) 


which  is  like  a cumulative  distribution  function  for  the  timescale. 
In  another  form  eq.  15  If 


6p  - 6w 


1.03  X 10“*^  e ~*0®58B2 

\ ‘ 


f e"*^/2 

X 


dz , 
(16) 


„ inM„  - m -0.78(?2 

8 

Here  can  be  related  to  This  equation  applies  only  when  the 

behavior  in  between  that  for  water  and  maximum  drag  reduction.  Also 
the  proportionality  factor;  should  be  a function  of  Reynolds 
number,  U*B/v.  Then  it  would  be  best  to  use  Ag  as  a characterization 
factor,  since  a ratio  would  eliminate  the  Reynolds  number  dependence. 


f 
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A 

8 


6-6  C <M  >0-70j.-.  0050B  ^ Jx  o“Z^/2 
_E _ P P E 

/ c“^^/2  ^jz 


(17) 


Is  a possible  ratio  where'  the  denominator  refers  to  maximuni  drag 
reduction.  Unfortunately  thin  requires  a standard  of  comparison 
for  each  U*D/v.  One  standard  which  can  be  deduced  from  previoua 
data  is  the  limit  when  is  negative  infinity. 

Consider 


«M-  «w 


0 . 7 8 g- . 0 8 5 0 6j7^ 


1 

■—  / 
^ \ 


(18) 


where  K is  the  proportionality  factor.  In  Patterson  and  Abernathy 
(1970),  <m2>/<M>2  was  found  to  be  appro: 'matcly  20  for  the  Polyox 
WSR  series.  This  gives  B ■ 1.73  and  sii.i  c <m2>/<M>,  the  weight 

V ‘ 

average  molecular  weight  is  8 x lO*'  we  can  also  get  <M>.  Then  for 
points  between  the  water  and  maximum  drag  reduction  cviivcs  we  can 
find  corresponding  to  -the  proper  fraction  of  6j^  - 6^.  For  such 
a wide  distribution  we  can  assume  that  the  effects  of  the  distribution 
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will  far  outwclf;h  any  Reynolds  number  dependence.  We  will  assume 
tbls  fraction  applies  to  the  flowrate  at  constant  Tor 

Patterson  and  Abernathy's  data 


Fraction  of 
(V«w> 

X 

«o 

Tl  . 

U*’/v 

0.01 

2.326 

5.7x10’ 

1.33xl0~^ 

2.25x10 

0.10 

1.282 

8.5x10^^ 

4. 5x10" 3 

C.6x102 

0.25 

.674 

2.97x10'" 

6.6x10“** 

4.5x103 

0.50 

0 

9.2x10^ 

8.3x10"3 

3.6x10*' 

0.75 

- .674 

2.85x105 

1.04x10"3 

2.94x10^' 

0.90 

-1.282 

1 X 10^ 

1.53x10“^ 

1.98x10® 

0.99 

-2.326 

1.49x10*' 

5.2x10“® 

5.7x10’ 

These  points  are  plotted  on  fig.  3 and  In  terms  of  Reynolds  number 
vs.  friction  factor  on  fig.  A.  The  concentr.it Ion  corresponding 
to  this  curve  is  approxim.itely  A ppm  by  weight  on  fig.  Similar 
curves  for  other  concentrations  can  he  prepared  by  multiplying  the 
fraction  of  (6j^  ~ 6^,)  in  the  table  by  C/4  where  C is  the  concentration 
in  parts  per  million.  Curves  for  various  multiples  of  the  rlnlmum 
curve  are  also  shown  on  fig.  3 and  4 and  compared  with  experimental  data 
on  fig.  4.  When  concentrations  are  much  less  or  much  more  than  the 
ifllnlmuin  to  reach  maximum  drag  reduction,  the  agreement  is  only  fair 
since  the  log  normal  distribution  is  only  an  approximation. 


-11- 


Aprll  6,  1973 
NSBrmac 


However,  the  recult  qualitatively  describes  the  behavior  quite  v;ell 

especially  if  the  straight  line  portiona  of  the  carves  on  fig.  3 

arc  extrapolated  to  the  water  and  maxitnun  drag  reduction  curves. 

We  see  that  many  observations  on  dilute  poljnner  solutions  are 

merely  duo  to  the  molecular  weight  distribution.  ! 

i 

To  interpret  tlie  Polyometcr  curves  we  see  that  the  onset 

point,  slope  of  the  curve  between  onset  and  maxlmuia  drag  reduction 

and  the  intersection  point  with  the  maximum  drag  reduction  line 

are  functions  of  concentration, average  molecular  weight  and  the 

parameter  B related  to  the  ratio  of  weight  average  to  number 

average  molecular  weight.  It  would  be  easier  to  relate  the  Intersection 

point  with  the  maximum  drag  line  with  tlie  appropriate  fresh  Polyox 

WSR  301  since  the  exact  nature  of  the  degradation  is  unknown. 

Another  possibility  is  to  use  the  appropriate  shear  rate  for  the 

\ 

test  and  relate  this  to  curves  like,  those  on  fig.  3.  As  we  can 
see  on  fig.  3 the  curves  depend  on  the  pipe  size. 

Fig.  5 and  6 chow  the  correlation  of  maximum  drag  reduction 
Intercept  with  concentration  of  fresh  Polyox  WSR  301.,  Using  fig. 6 ' 

for  the  A8-lnch  tunnel  test  we  obtain  the  following  for  a speed 
of  8.5  meters  per  second 
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Test 

Equivalent 

Fresh 

Concentration 

7627 

> 10 

7629 

10 

7632 

7 

7633 

5 

7635 

3 

7639 

1.8 

7641 

2.5 

This  equlvalonC  would  be  different  for  each  speed  since  the 
onset  point  is  not  equivalent  to  the  original  Polyox.  For  example 
text  7633  corresponds  to  an  onset  shear  rate  of  4 x 10^  scc“^. 

On  fig.  5 the  onset  point  is  plotted  for  this  range  and  at 
4 X 10^  sec  we  obtain  1.5  ppm  compared  to  5 ppm  from  the  comparison 
at  a shear  rate  of  8 x 10**  sec~*.  Obviously  the  concentration 
is  greater  than  any  of  these  based  on  fresh  Polyox  WSR  301  but  J 
the  molecular  weight  average  is  less.  We  can  estimate  the  j 

average  molecular  weight  and  concentration  if  vv  assume  B is 
unchanged.  This  is  reasonable  only  because  B is  quite  large  anyway 
and  the  probable  increase  in  the  weight  average  to  number  average 
molecular  treight  ratio  would  not  change  3 by  much.  Let  us  further 
assume  that  the  concentration  remains  constant  so  only  molecular 
weight  changes  and  plot  the  results  on  fig.  7.  The  minimum 


I 
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concentration  to  reach  maxlmuni  drag  reduction  for  various  weight 
average  molecular  weights  Is  found  to  be  as  follows: 

Mw 

8 X 10^ 

A X 10^ 

2 X 10^ 

1 X 10® 

5 X 10® 

The  log  normal  model  predicts  minimum  concentrations  to  reach 

maximum  drag  reduction  in  agreement  with  Patterson  and  Abernathy. 

In  addition  the  apparent  "onset"  shear  rate  increases  with  decreasing 

\ • 

average  molecular  weight  and  decreases  with  concentration  increase 
also  in  agreement  id.th  their  data.  We  also  find  the  same  behavior 
in  stretching  out  the  curves  when  the  "onset"  Is  within  the  laminar 
flow  regime  for  the  pipe  used  as  shown  on  fig.  3.  Much  of  the 
"unusual"  behavior  of  the  polymer  solutions  are  therefore  only  a 
rcrfult  of  the  molecular  weight  distribution.  Fig.  7 can  be  used 
to  estimate  the  weight  average  molecular  weight  for  the  48-lnch 
tunnel  test. 


Mlnlrauia  Concentration 
for  Maximum  Drag 
Reduction 

A ppm 

6.9  ppra 

11.7  ppm 

20.3  ppm 

3A . 7 ppm 
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The  approximate  results  are  as  follows: 


Test 

M 

w 

7627 

8 X 

10^ 

7629 

3 X 

10^ 

7632 

2 X 

10^ 

7633 

1.5 

X 10*' 

7635 

1 X 

10^ 

7639 

0.7 

X 10^ 

7641 

0.9 

X 

O 

Effects  other  than  drag  reduction  are  found  for  Increased 
concentrations  for  fresh  solutions  and  arc  not  Included  In  this 
analysis. 

\ » 

We  can  see  from  this  analysis  that  a polymer  ocean  test  can 
be  characterized  by  the  equivalent  of  fresh  polymer  or  the 
average  molecular  weight.  Either  of  these  can  be  estimated  as 
described  herein.  In  a separate  attaclunent  to  this  report  a 
"cookbook"  method  for  preparing  these  estimates  is  given. 
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Preparation  of  Polyorreter  CurvoB  from  Lof;  liomal  IJodel 

We.  assume  that  for  Polyox  WSU  301  the  vcl p,ht  avernge 
molecular  vclght  Is  8 x 10^,  the  ratio  of  weight  average  to 
number  average  molecular  weight  Is  20  and  the  minimum 
concentration  to  reach  maximum  drag  reduction  line  Is  4 ppm. 
Polntn  on  the  master  curve  are  an  follows: 


fraction 
•»  fl 

X 

Shear-rate, 
sec“^  “ U*^/v 

.99 

-2.326 

5,7  X 10^ 

.90 

-1.282 

1.98  X 10^ 

.75 

- .674 

2.94  X 10^ 

.50 

0 

3.6  X lo'* 

.333 

+ .432 
» ' 

9.5  X 10^ 

.25 

.674 

4.5  X 10^ 

.20 

.842 

2.73  X 10^ 

.15 

1.036 

1.5  X 10^ 

.10 

1.282 

6.9  X 10^  . 

.08 

1.405 

4.6  X 10^ 

.06 

1.555 

3.05  X lO^’ 

.0^ 

1.751  ' 

1.63  X 10^ 

.02 

2.054 

6.4  X 10 

.01 

2.126 

2.75  X 10 

f 


r 


j 

I 

I 


I 

1 

* 


I 
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Case  I - No  degradation  only  concentration  is  changed. 

a.  Multiply  the  numbers  In  the  fraction  column  by  the 
concentration  In  p])m  divided  by  U,  Disregard  all 
numbers  greater  than  1. 

b.  Plot  the  resu-ting  fraction  and  corresponding  shear  rate. 

At  the  constant  shear  rate,  read  the  flow  rates  for 
water,  Q^,  and  maximum  drag  reduction,  Q^.  Calculate 

Qp  “ Q„  + fi  (Qp  - Qw)- 

Plot  Qp  at  the  shear  rate. 

c.  When  sufficient  points  have  been  located,  draw  a straight 
line  through  the  points  in  the  central  region  between  the 
water  and  maximum  drag  reduction  curves  and  extrapolate 
the  line  to  these  two  curves.  Disregard  the  points  not  on 
this  line. 

Case  II  - Typical  polymer  ocean  test , concentration  is  constant  but 
molecular  weight  is  changed. 

a.  Calculate  the  ratio  8 x 10^  to  the  weight  average  molecular 
weight  of  the  curve  to  bo  plotted.  Calculate  tlie  following: 

S H (8xl0^/li^,)°*2® 

W B 

b.  The  new  minimum  concentration  to  reach  maximum  drag  rcductic.n 
is  AS.  The  fraction  numbers  are  multlj)lled  by  concentration 
divided  by  AS  and  the  shear  rate  numbers  are  multiplied  by 

W.  Form  a new  table. 

c.  I’roceed  r.s  in  case  I,  b and  c. 
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